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Восстановление томографической 
информации в системах визуализации гамма- 
излучения с кодированными апертурами 


В статье предложен метод восстановления томографической информации по результатам единственного 
наблюдения трехмерного распределения источников системой визуализации гамма-излучения с 
кодированной апертурой. Приведены результаты моделирования поля счетов и восстановления 
распределений источников в Плоскостях на различном удалении. Установлено соотношение ДлЯ 
расстояний между плоскостями разрешения. 


Метод визуализации гамма-излучения с помощью кодированных апертур представ- 
ляет собой двухступенчатую процедуру пространственного кодирования положения 
источника в поле зрения. На первой стадии направление падающего гамма-излу- 
чения кодируется с помощью расположенных перед позиционно-чувствительным 
детектором прозрачных и поглощающих элементов, на второй стадии зарегистрирован- 
ная позиционно-чувствительным детектором картина (так называемая «тенеграмма») 
подвергается декодированию для получения изображения поля зрения. 

Приборы для визуализации гамма-излучения, созданные на основе этой концепции 
пространственного кодирования направления прихода гамма-лучей, называют «сис- 
темами с кодирующими масками» или «системами с кодированными апертурами» — 
СКА. Характерными особенностями конструкции этих систем является наличие 
кодирующей маски, состоящей из прозрачных и поглощающих гамма-излучение 
элементов, расположение которых в апертуре подчиняется строго определенным 
математическим закономерностям, и наличие позиционно-чувствительного дете- 
ктора (ПЧД), пространственное разрешение которого должно быть согласовано с 
размером элемента кодирующей маски. 

Фотоны источника гамма-излучения в поле зрения (ПЗ) проектируют тень 
маски на детектор. В течение единственного наблюдения детектор собирает «тене- 
граммы» от каждого источника в ПЗ, при этом информация о позиции источника 
кодируется сдвигом тенеграммы, а интенсивность источника — мощностью тенеграммы. 
Конфигурация кодирующей маски разрабатывается так, что тень маски, проекти- 
руемая источником в любой позиции в поле зрения уникальна и, следовательно, 
записанная тенеграмма может быть использована для восстановления положения 
источника. Принципы визуализации гамма-излучения с помощью кодированных апертур 
и некоторые способы конструирования кодирующих масок достаточно подробно 
изложены в [1], [2]. 

Счетная информация, зарегистрированная позиционно-чувствительной плоскостью 
в течение времени наблюдения, представляет собой суперпозицию тенеграмм от 
всех источников в поле зрения. При этом тенеграммы от источников, находящихся 
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на различных расстояниях от ПЧД и кодированной апертуры, имеют различный 
масштаб в соответствии с законами геометрической оптики [3]. Принцип регистрации 
тенеграмм от источников на различных расстояниях от ПЧД показан на рис. 1. 
Фокальной плоскостью в данном случае называется плоскость, источники которой 
проецируют базовый фрагмент маски на всю чувствительную поверхность ПЧД. 

Таким образом, трехмерная реконструкция распределения источников может 
быть получена путем исследования данных детектирующей плоскости для теней, 
соответствующих источникам на различных расстояниях. Отметим, что подробный 
анализ тенеграмм возможен только в случае применения в качестве позиционно- 
чувствительной плоскости детектора непрерывного типа (гамма-камера) или пиксе- 
лированного детектора, поперечный размер индивидуальных кристаллов которого 
значительно меньше размера тени от элемента маски. 


Кодированная 
апертура 


Источник вне фо- 
кальной плоскости 


Источник в фо- 
капьной плоскости 


Рисунок 1. — Принцип регистрации томографической информации в СКА 


Цель работы -— разработка метода восстановления томографической информации 
по результатам однократного наблюдения объемного источника гамма-излучения 
системой визуализации гамма-излучения с кодированной апертурой. 

При прохождении потока гамма-квантов, излучаемых источниками в поле зрения, 
через апертуру системы визуализации происходит его пространственная модуляция 
вследствие взаимодействия гамма-квантов с материалом поглощающих элементов 
кодирующей маски. В результате на ПЧД падает поток, отображающий некоторую 
часть маски. При регистрации потока излучения детектором события взаимодействия 
записываются в некоторый дискретный массив счетов, отображающий координаты 
точки взаимодействия в плоскости ИЧД и выделенную при взаимодействии энергию. 
Этот массив мы называем полем счетов ПЧД. Количество актов взаимодействия 
гамма-квантов с кристаллом ПЧД за время наблюдения пропорционально плотности 
потока, поэтому для получения картины тени (тенеграммы) необходимо провести 
интегрирование количества актов взаимодействия в элементах ПЧД, являющихся 
проекциями элементов кодирующей маски. 

На рис. 2 показан пример поля счетов ПЧД при наблюдении трех протяженных 
источников, находящихся в плоскостях на различных расстояниях от кодирующей 
маски. Для модели использована кодирующая маска размерностью МхМ = 17х17 
элементов МОКА [2]. 

В случае использования в качестве ПЧД системы визуализации детектора 
непрерывного типа (гамма-камеры) размеры элемента интегрирования определяются 
по законам геометрической оптики [3]: 


ана [1 =а-, (1) 
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где 4 — линейный размер элемента интегрирования ПЧД; а — размер элемента 
маски; /— расстояние ПЧД-маска; Ь — расстояние маска — плоскость источников; т — 


коэффициент увеличения. 
тра М | 


та М 


тра М 


Ш М 


Рисунок 2 — Пример поля счетов ПЧД при наличии источников 
в поле зрения на различных расстояниях 


Для восстановления распределения источников излучения в поле зрения системы 
визуализации чаще всего используется метод корреляции [1]. В этом случае про- 
странственное распределение источников О' получают в соответствии с формулой: 


№М-—1М-1 


О, № ЕЁ, "бк › (2) 


где С — обрабатывающая функция, Р — матрица счетов, М, № — размерность 
кодирующей маски. Матрица счетов Р является результатом интегрирования поля 
счетов по площади тени элемента маски. Размерность матрицы Р совпадает с 
размерностью базового фрагмента маски. 

Источник излучения, находящийся в фокальной плоскости, создает на поверх- 
ности детектора тенеграмму в виде фрагмента маски размерностью, совпадающей с 
размерностью базового фрагмента маски. Источники в удаленных плоскостях про- 
ецируют больше элементов, чем размерность базового фрагмента маски. Для полу- 
чения матрицы счетов, обеспечиваемой источниками в удаленных плоскостях, поле 
счетов ПЧД необходимо интегрировать только по некоторой части, по площадкам, 
размеры которых совпадают с размером тени элемента. На рис. 2 показаны размеры 
площадей интегрирования для получения матрицы счетов, обеспечиваемых источни- 
ками в различных плоскостях. Коэффициент увеличения то соответствует фокальной 
плоскости. Соотношение коэффициентов увеличения: то > т, > т2. Результат 
восстановления источников в трех плоскостях при различных значениях коэффициента 
увеличения показан на рис. 3. 

Источники в плоскостях ближе к ПЧД, чем фокальная плоскость, создают на 
поверхности ПЧД тенеграмму с количеством элементов, меньшим, чем размерность 
базового фрагмента. Таким образом, ПЧД регистрирует недостаточно информации 
для восстановления распределения источников. Поэтому распределение источников 
может быть восстановлено только в небольшом диапазоне удалений от фокальной 
плоскости. 
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в) 
Рисунок 3 — Результат восстановления источников в различных плоскостях: 
а) фокальная плоскость 770; 6) т! < то; в) т2 < т! 


В результате суперпозиции тенеграмм от источников в различных плоскостях 
на восстановленных изображениях просматриваются также и источники в соседних 
плоскостях, однако их уровень ниже уровня источников в восстанавливаемой 
плоскости. На рис. 3 а) кроме источника в виде ромба внизу слева отчетливо виден 
центральный крест из соседней плоскости. На рис. 3 6) отчетливо просматриваются 
следы источника внизу слева, находящегося в фокальной плоскости. Влияние соседних 
плоскостей особенно значимо для источников в наиболее удаленных плоскостях. 
Так, на рис. 3 в) источник в виде ромба вверху справа восстановлен только в виде 
отдельных точек. Очевидно, что такие результаты восстановления требуют дополни- 
тельной обработки для уточнения картины распределения источников в восстанавли- 
ваемой плоскости. 

Для очистки изображения источников от влияния соседних плоскостей можно 
применить методы субтрактивной деконволюции [4]. Проведение такой обработки 
требует знания функции рассеяния точки (ФРТ) по глубине, т.е. зависимость отклика 
системы визуализации от расстояния между источником и кодированной апертурой. 
Исследованию таких зависимостей посвящены наши предыдущие работы [3]. 
Другим методом, успешно применяемым для решения задачи очистки от влияния 
помех и шумов, является метод максимума энтропии [5]. 

Представляет интерес величина томографического разрешения, обеспечиваемого 
методом интегрирования поля счетов. Естественно считать, что результаты восста- 
новления источников могут различаться, если различаются результаты интегрирования 
поля счетов. Различие результатов интегрирования обеспечено различием площадей 
элементов интегрирования (тени элемента). Если поле счетов ПЧД записывается с 
дискретностью 9, то размеры элементов интегрирования могут минимально разли- 
чаться на эту величину. Используя соотношение (1), получим 


а:(т-т)=б, (3) 


где а — размер элемента маски, то — коэффициент увеличения для фокальной 
плоскости, т — коэффициент увеличения для плоскости восстановления распределения 
источников. 

Используя определение коэффициента увеличения, из (3) получим соотношение 
для расстояний между плоскостями разрешения: 


ры. 2 

ь Ба: 
Здесь Бо — расстояние от кодирующей маски до фокальной плоскости, В — 
расстояние от кодирующей маски до плоскости восстановления, /-— расстояние между 
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маской и ПЧД. Очевидно, что геометрические параметры системы визуализации 
гамма-излучения имеют определяющее влияние на ее томографические возможности. 
Предложенный в статье метод восстановления томографической информации в 
системах визуализации гамма-излучения с кодированными апертурами путем интегри- 
рования поля счетов по площади тени элемента маски дает возможность восстановить 
распределение источников излучения в плоскостях, находящихся на различном удалении 
от ПЧД системы. Особенность и ценность этого метода состоит в том, что всю информа- 
цию для восстановления томографических изображений система визуализации получает в 
течение единственного наблюдения трехмерного источника с одного ракурса. Получа- 
емые в результате восстановления изображения оказываются зашумленными в результате 
влияния соседних слоев. Однако методы очистки изображений позволяют устранить 
шумы и обеспечить точную картину распределения источников. Приведенное в статье 
соотношение для расстояний между плоскостями разрешения позволяет установить 
геометрические размеры источников в плоскостях, что необходимо для восстановления 
трехмерного распределения источников по двумерным томографическим срезам. 
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В.Ю. Плахотвак 

Вищновлення томограф!чно! 1нформаци в системах взуалвацй гамма-випромнювання 

з кодованими апертурами 

У статт!г запропонований метод вдновлення томографчно! 1нформацй за результатами единого 
спостереження тривимфного розподллу джерел системою взуалзацй гамма-випроминювання з кодованою 
апертурою. Наведен! результати моделювання поля л!чби 1 вдновлення розподшив джерел в площинах на 
разних вдстанях. Встановлено сшввдношення для вадстаней мж площинами розрзнення. 
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